o 



THE BRIGHTER SIDE OF SUB-MM SOURCE COUNTS: A 
SCUBA SCAN-MAP OF THE HUBBLE DEEP FIELD 

COLIN BORYS^ SCOTT CHAPMAN*, MARK HALPERN^DOUGLAS 

SCOTT* 

^Department of Physics & Astronomy, University of British Columbia, Vancouver 

BC CANADA V6T IZl 
^^ ' E-mail: borys@physics.ubc.ca 

^'- 

^\\ ■ * Observatories of the Carnegie Institution of Washington, Pasadena CA USA 

^. 91101 

j-^ ■ We present an 11 x 11 arcminute map centred on the Hubble Deep Field taken 

^ ■ at 850 /im with the SCUBA camera on the JCMT. The map has an average one- 

^y\ ' sigma sensitivity to point sources of about 2.3 mjy and thus probes the brighter 

end of the sub-mm source counts. We find 7 sources with a fiux greater than 

^_^ , 9mJy (~ 4cr), and therefore estimate A^(> 9mJy) = 208 _72 degree ~-'^. This result 

^ I is consistent with work from other groups, but improves the statistics at the bright 

-^ ■ end, and is suggestive of a steepening of the counts. 

^: 

^H ■ 

^\ • 1 Introduction 

o: 

^D ■ Observations using the Submillimetre Common User Bolometer Array 

^^ I (SCUBAtl) have been revolutionizing our understanding of the importance 

^ r^ . of dust in galaxies at high redshift. Already several very deep integrations 

Qh' have been carried out on single SCUBA fields down to the confusion limit. In 

Q . order to learn more about source counts (and hence to constrain models), the 

?-H I next step is to search for brighter objects over somewhat larger fields. 

C/) . The population oLbright sub-mm sources is currently not well under- 

f^ ' stood. Current modelsHooB for source counts in the sub-mm have been able 

^ . to account for the observed sources by invoking evolution which follows the 

(1 + z)^ form required to account for IRAS galaxies at 60 /im, and the powerful 
radio-galaxies and quasarsQ. Euclidean models with no evolution have a slope 
C^ I of roughly —1.5 {Ng oc Su ), which cannot possibly account for the sources 

observed. With reasonable evolution (in IRAS-motivated models), the counts 
steepen sharply at the lO's of mJy sensitivity level to roughly S~'^-^. At the 
source detection level of previous work, typically 5 mJy, there is little variation 
between the models, and we are well into the steep counts regime. _ 

However, at the 10-30 mJy sensitivity level various evolutionary modelaO 
show more parameter dependence. Given a possibly wide range of galaxy 
types contributing to the source counts at the bright end, the actual counts 
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may deviate from a simple parametric model dramatically. Furthermore, cos- 
mological, as well as evolutionary parameters, play a role here. 

2 Data collection and analysis 

The standard scan-map observing strategy is to use multiple chop throws in 
two fixed directions on the sky so that an Emerson deconvolution technique 
can be applied. This may be appropriate for galactic plane mapping, where 
structure appears on all scales, but we seek a simpler approach which is opti- 
mized for finding point sources. The disadvantage of the standard approach is 
that the off-beam pattern gets diluted (making it more challenging to isolate 
faint sources), and the map noise properties are more difficult to understand. 
Our approach was to use a single chop direction and throw fixed on the sky. 
The result is a map that has, for each source, a positive and negative signal. 

A total of 61 scans were obtained, 31 in early 1998 and 30 more a year 
later. The data were first analysed using SURF, but because the data were 
taken in a non-standard mode, we found it necessary to write custom software. 
In this way, we were able to isolate and remove large scale features in the maps, 
and estimate the per pixel noise level by a careful accounting of the per-bolo 
noise and the frequency at which it sampled a particular pixel. 

Sources were found by using a model for the dual beam pattern and fitting 
it, in a least-squares sense, to each pixel in the map. A total of 7 sources were 
found that had a peak to error-of-fit ratio greater than 4. We also rotated 
the dual beam pattern by ninety degrees and used that as a model; Only 2 
sources were found, one of which is associated with the brightest source in 
the map. Monte-carlo simulations suggest that the other false positive is not 
unexpected. As an additional check, the 1998 data were compared with the 
1999 data to ensure that each source was evident in both halves of the data. 

The calibration was determined by fitting the beam model to observations 
of standard calibrators taken during the run. We checked our procedure by re- 
analysing the BargerB et al. HDF/radio fields, which are a set of jiggle maps 
within the HDF fianking fields. In fact this led to a correction in our original 
fiux estimates. Figure 1 presents our final map, and an updated source count 
plot is given in Figure 2. 

3 Conclusions and followup work 

We have significantly improved statistics at the ~ 10 mJy level, and suggest 
there is some indication of a steepening of the counts there. Because these 
sources are brighter than typical SCUBA detections, they will be relatively 
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Figure 1. The 850 /im HDF map. The map size is roughly 12 x 12 arcminutes, although only 
the central 11 x 11 was used in order to avoid the noisy edges. The chop is 40 arcseconds 
and is roughly east- west. The circles outline the 7 sources detected in our survey and the 
dark crosses highlight sources found in the Barger et al. (2000) and Hughes et al. (1998) 
work. Since these surveys went deeper, not all can be detected in the scan map. They also 
cover a smaller area, and hence cannot see several of our sources. 



easy to follow up at other wavelengths. Preliminary analysis already indi- 
cates good correlation with jaJy radio sources, but little indication of optical 
counterparts (as found in other studies). 

Although not discussed here, we are also investigating sources in the 
450 /im scanmap. However, this is more arduous given the lower sensitiv- 
ity of SCUBA with this filter. Finally, a similar analysis is being conducted 
on data taken of 50 square arminutes of the Groth strip/Westphal region. 
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Figure 2. The 850 /xm source counts. The crosses are data from Blain et al. (1999) and the 
open circles that from a recent cluster survey by Chapman et al. (2000) The dashed line is 
a simple two power-law fit to the previous counts. Our new estimate, with Bayesian 68% 
confidence limits, is given by the solid circle. 
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